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Abstract 

We derive van der Waals-London and Casimir forces by calculat- 
ing the eigenmodes of the electromagnetic field interacting with two 
semi-infinite bodies (two halves of space) with parallel surfaces sep- 
arated by distance d. We adopt simple models for metals and di- 
electrics, well-known in the elementary theory of dispersion. In the 
non-retarded (Coulomb) limit we get a d _3 -force (van der Waals- 
London force), arising from the zero-point energy (vacuum fluctua- 
tions) of the surface plasmon modes. When retardation is included 
we obtain a d~ 4 - (Casimir) force, arising from the zero-point energy 
of the surface plasmon-polariton modes (evanescent modes) for met- 
als, and from propagating (polaritonic) modes for identical dielectrics. 
The same Casimir force is also obtained for "fixed surfaces" boundary 
conditions, irrespective of the pair of bodies. The approach is based 
on the equation of motion of the polarization and the electromagnetic 
potentials, which lead to coupled integral equations. These equations 
are solved, and their relevant eigenfrequencies branches are identified. 

Key words: Electromagnetic eigenmodes in matter; van der Waals-London 
and Casimir forces; surface plasmons and plasmon-polaritons 
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1 Introduction 

The Casimir force was originally derived by estimating the zero-point energy 
(vacuum fluctuations) of the electromagnetic field comprised in-between two 
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ideal, perfectly reflecting, semi-infinite metals (two halves of space) separated 
by distance d[T] As it is well-known, it goes like <i~ 4 for distances greater than 
the characteristic electromagnetic wavelengths of the bodies (plasmon "wave- 
lengths"). Further on, the calculations have been cast in a different form, by 
resorting to the fluctuations theory, [21 13] and a d _3 -force has been obtained 
for the non-retarded (Coulomb) interaction, which corresponds to the van 
der Waals-London force. The matter polarization is usually represented in 
this case by a dielectric function. Recently, there is a renewed interest in this 
subject, motivated, on one hand, by the role played by plasmons, polaritons 
and other surface effects arising from the interaction between the electromag- 
netic field and matter and, on the other hand, by the querries related to the 
applicability of a dielectric function for discontinuous bodies. [1]-[2D] We re- 
port here on a different investigation of these forces, based on the calculation 
of the eigenfrequencies of the electromagnetic field interacting with matter. 

We assume a simple model of matter, consisting of mobile particles with 
charge — e and mass m, moving in a rigid neutralizing background, and sub- 
jected to certain forces. Such a model is reminiscent of the well-known jellium 
model of electron plasma, though it is generalized here to some extent. In 
the presence of the electromagnetic field matter polarizes. We leave aside 
the magnetization (we consider only non-magnetic matter) and relativistic 
effects. We represent the small disturbance in the density of the mobile 
charges as 5n = —ndivu, where n is the (constant) concentration of the 
charges and u is a displacement field in the positions of these charges. The 
charge disturbance is therefore p = endivu. This representation is valid for 
Ku(K) <C 1, where K is the wavevector and u(K) is the Fourier transform 
of the displacement field. 

For homogeneous and isotropic matter the displacement field obeys an equa- 
tion of motion which can be taken of the form 

mii = — eE — eE — mu^u — mju , (1) 

where E is the (internal) electric field, E is an external electric field, u is a 
frequency parameter corresponding to an elastic force and 7 is a dissipation 
parameter. Making use of the temporal Fourier transform we get 

uH = -- \- — (E + E ) (2) 

(where we dropped out the argument uj of the electric fields). On the other 
hand, from Maxwell's equation divE = inendivu, we get the (internal) elec- 
tric field E = 4imeu (equal to — 47rP, where P is the polarization). Making 
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use of equation (J2J) we get the dielectric function 

eH = 1 - 2 j— (3) 

from its definition E = e(E + E ), where u pi given by u 2 = 47me 2 /m, is 
the plasma frequency. The dielectric function given by equation (J3|) is well 
known in the elementary theory of dispersion. [21] It proves to be a fairly 
adequate representation for matter polarization in various bodies. We can 
view u p ,u}q and 7 as free parameters, thus being able to simulate various 
models of matter. For uiq = 7 = we get the well-known dielectric function 
of an ideal plasma; if ujq = we have the dielectric function of the optical 
properties of simple metals for uj 7 (Drude model), and the dielectric 
function corresponding to the static (or quasi-static) currents in metals for 
uj <C 7; for ujq ^> Up we have a dielectric function of dielectrics with loss; and 
so on. 

In addition, making use of equation (J2J), we can compute also the electric 
conductivity a, from its definition j = cr(E + E ), where j = — enii is the 
current density. We get the well-known conductivity 

a(uj) = -p 1 — 2— — , (4) 

47T U 2 - UJ 2 + lOJ^f 

whence, for instance, the static conductivity for metals a = Up/in'j; param- 
eter 7 can be viewed as the reciprocal of a damping time r (or relaxation 
time, or lifetime), 7 = 1/r, and we get the well-known static conductivity 

o = ne 2 r/m'~f. 

Therefore, the equation of motion (0Q) turns out to be an adequate starting 
point for representing the matter polarization. However, we must note that 
for dielectrics, which may imply oscillations in localized atoms (in our model 
through the frequency ojq), the classical dynamics assumed here turns out 
to be inadequate in the retarded regime, and a quantum treatment is then 
required. 

In the non-retarded limit the electric field E in equation (0Q) is given by the 
Coulomb law, i.e. E = —grad&, where $ is the static Coulomb potential 
arising in matter. The latter depends on the charge disturbance p = —eSn, 
therefore on u. Then, it is easy to see that the equation of motion (DQ) leads 
to an integral equation for the displacement field u. Its eigenvalues give the 
plasmon modes. For retarded interaction, the electric field E in equation 
(PQ) is given by the vector potential A and the scalar potential $ trough 
E = — grad&. Making use of the radiation (Kirchhoff) formulae, these 
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potentials can be expressed as integrals containing the displacement field 
u (through the charge and current densities), and we get again an integral 
equation for u. Its eigenvalues give polariton-like modes. The use of integral 
equations in treating the electromagnetic field interacting with matter was 
previously indicated in connection with the so-called Ewald-Oseen extinction 
theorem. [22] We have applied this approach to a semi-infinite (half-space) 
body, as well as to a slab of finite thickness. [23] In this case, beside the bulk 
displacement field, there appears a surface displacement field also, and the 
integral equations couple these degrees of freedom. We have solved these 
coupled integral equations and computed bulk and surface plasmons and 
polaritons, dielectric response, reflected, refracted and transmitted fields, and 
derived generalized Fresnel relations. We employ the same procedure here 
for two semi-infinite bodies (two halves of the space) separated by distance 
d, in order to get the electromagnetic eigenfrequencies and to derive van der 
Waals-London and Casimir forces. We do it in two steps: first, for static 
Coulomb (non-retarded) interaction (valid for wavelengths much longer than 
the characteristic size of the bodies) and, second, for retarded interaction. 



2 Surface plasmons. van der Waals-London forces 

We consider two semi-infinite bodies (two halves of space) with parallel sur- 
faces in the (x, y)-plane, separated by distance d. The bodies occupy the 
regions z < —d/2 and, respectively, z > d/2. We take two displacement 
fields u 12 , giving rise to two charge disturbances Sn 12 = — ni ^div u lj2 . We 
consider first the equation of motion for an ideal plasma. In general, we leave 
aside the dissipation (parameter 7 in equation (pQ)), which is irrelevant for 
our discussion. The equation of motion reads 



miii = grad J dH U ( 



R- R 



nidivuiCR ) + n 2 <i2t>u 2 (R 



(5) 



and a similar equation for u 2 , which can be obtained from equation ([5]) by 
interchanging the labels 1 and 2 (1 < — > 2); U(R) = e 2 /R in equation ((SJ) 
is the Coulomb interaction. Since we are interested in the eigenmodes, we 
leave aside the external field E . We use R = (r, z) for the position vector 
R, where r = (x,y), and the representation 

ui, 2 =(vi >3 ,wi, 2 )0(±3-d/2) (6) 

for the displacement fields, where 9(z) = 1 for z > and 9(z) = for z < 
is the step function; the ± sign is associated with labels 1 and 2, respectively. 
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The divergence in equation (EJ) can now be written as 

divv 1}2 + -j^ J 6(±z - d/2) + w 1>2 {±d/2)6(±z - d/2) , (7) 

where wi t2 (±d/2) means 101,2(1", 2 = ±d/2). We notice in equation the 
(de) polarization charge arising at the surfaces z = ±d/2. We employ Fourier 
representations of the form 



vi, 2 (r,z;t) = / da;vi, 2 (k,. 

1, J 



z; w)e ikr e^ (8) 



and similar ones for W12, and use the Fourier transform 



Vr 2 + z 2 ^ k 

for the Coulomb potential. Then, we notice that equation j5|) implies that 
vi )2 are parallel with the wavevector k (in-plane "longitudinal" modes), and 
zku*!^ = We use this latter relation to eliminate wi t 2 from the equations 

of motion. In addition, we introduce the notation v\ j2 = kvi^/k. Then, it is 
easy to see that equation (JHJ) yields two coupled integral equations 

^2 - f°° dz e- k \ z ~ z '\ih + ^ f°° dz -9-e~ k \ z - z '\^+ 
u vi - 2 j d/2 tt2;e u i + 2fc Jd/2 ttZ a 2 ,e a/^ 

: jk r-d/2 , / k r z _ z /\ u>l r—d/2 



(10) 



^ = * Q dz'e*(>-*'\ - f J- dzW-*%+ 

+ &^tie-*Mv 2 + £^dz&e->^*\%, z<-d/2 , 

where oj\ 2 = 4nni t2 e 2 / m and we dropped out the arguments u>, k. Integrating 
by parts in equations (fTOl we obtain a system of two algebraic equations 

(to 2 - ul)v x = -\e~ kz [oj 2 e kd ' 2 Vl {d/2) - uj 2 e- kd ' 2 v 2 {-d/2)] , z > d/2 , 

(uj 2 - uo 2 )v 2 = \e kz [uj\e~ kd l 2 v x {d/2) - uj 2 e kd / 2 v 2 (-d/2)] , z < -d/2 . 

(11) 

We can see that in this non-retarded limit the two bodies are coupled only 
through their surfaces. 
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For Vi(d/2) = v 2 (— d/2) = in equations (TTTT) we get the bulk plasmons 
uj = Ui ; 2- Making z = ±d/2 in equations ffTTl) we get the system of equations 
for the surface modes. The corresponding dispersion equation is given by 

- 2 " ^l) (- 2 " - = • (12) 

For d = we obtain the surface plasmon of a metallic interface given by 



UJ 
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2 - 1 'u)\ + cjf), while for d — > oo we get the surface plasmons uj = uJ\ >2 /\/2 
for free (uncoupled) surfaces. If the body labelled by 2 for instance is a 
dielectric, then uj 2 in the second equation IjHl) is replaced by uj 2 — uj 2 . In the 
limit ujo ^> UJ2 and for d = we get the surface plasmon uj = uji/y/1 + e 2 , 
corresponding to a dielectric-metal interface, where £2 = 1 + k>IA t 'o- For two 
identical metals c^i = uj 2 = uj p we get the surface plasmons given by 

uj 2 = \uj 2 v (1 ± e~ kd ) . (13) 

They are identical with the surface plasmons of a plasma slab of thickness d. 
These are well-known results. [24|-|31j 

Let us label by a all the eigenvalues Q a of the system of equations (fTTl) . We 
compute the force acting between the two bodies by 

d 



3d 

a 

where we recognize the zero-point energy of harmonic oscillators. Although it 
can be included straightforwardly, it is easy to see that the temperature plays 
no significant role, so we may neglect the temperature effects, as usually. We 
may also leave aside the bulk plasmons, since they do not depend on the 
distance d. We are left with the two surface modes Jli i2 given by equation 
(112]) , labeled by wavevector k. We can see that these eigenvalues are function 
of kd, so the force depends on distance d as F ~ 1/d 3 . As it is well-known, 
such a force between two bodies implies an inter-atomic interaction ~ 1/R 6 , 
where R is the distance between two atoms. This is the well-known van der 
Waals-London interaction.! 321 



We compute here the force F for the eigenvalues given by equation (TT3l) , i e. 
for two identical plasmas (metals). Equation (TTill gives a force 



hujp 



dx ■ x 2 e~ x - - (15) 



8irV2d?J \Vl-e- x Vl 

per unit area. The integral in equation (1151) is ~ 4, so we get F ~ huj p /2ir\^d 3 . 
In like manner we can compute the force between two (identical) dielectrics, 
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by replacing u 2 in equation ([TBI by u 2 — u 2 and taking the limit u> ^> u p . 
The result is a much weaker force F = hWp/128ujQd 3 . It can also be written 
as F = huj (e — l) 2 /128<i 3 , where e ~ 1 + uj 2 /uq is the (static) dielectric 
function in the limit u> <C ojq. The same result is obtained by making use 
of the formulae given in Ref. [32] for non-retarded interaction within the 
framework of the fluctuations theory (equation 82.3 p. 343 in Ref. [32]). 
Making use of the eigenvalues given by the roots of the dispersion equation 
(fl~2l) . we can compute in the same manner the force acting between two dis- 
tinct bodies. For instance, we can consider a dielectric-metal pair and get 
straightforwardly the force F = fauJitu%/327r^/2LU 2 l d 3 , where u)\ belongs to the 
metal and a>2, ooq represent the dielectric. 



3 Surface plasmon-polariton modes. Casimir 
force 

We pass now to the retarded interaction. The electric field in equation (pQ) is 
given by E = — — grad& } where A is the vector potential and $ is the 
scalar potential. These potentials are given by 

A(r, z;t) = lf dr> j d ^;t-R/c) (lg) 



and 



S(r,z;t) = / dv' j dz' p{r '> Z ' ]t R R,C) , (17) 



where 

j = -eni(vi, w^z - d/2) - en 2 (v 2 , w 2 )6(-z - d/2) (18) 
is the current density, 



p = eni (diuvi + 9(z - d/2) + Wl (d/2)5(z - d/2)+ 
+en 2 (divv 2 + 6(-z - d/2) + w 2 (-d/2)5(z + d/2) 



(19) 



is the charge density and R = y/fr — r') 2 + ( z ~ z ') 2 - We use the Fourier 
representations given by equation (JSj) and the Fourier transform|33j 



^ 2m 



Vr 2 + z 2 ^ k 
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J^—e^e^ , (20) 



where k = y ^ — k 2 . Then we compute the electric field from the potentials 
given by equations (fT6l) and (fTTj) and use equation (pQ) for oj$ = 0, 7 = 0, E = 
in order to get integral equations for vi 2 , ti>i,2- We define the wavevector 
k± of magnitude k and perpendicular to the wavevevctor k, and introduce 
the notations t>i 2 = kv 12 /A;, v± 2 — kj_v lj2 /fc. Doing so, we get the first set 
of integral equations 



V 



'I ~~ 2c 2 k Jd/2 



lj-dz'e-^%i{z') -§-j:^ dz > e i«\*-A v ^ z i) , z < _rf/ 2 



u 2 — 2c 2 k Jd/2 ' 

(21) 

where we dropped out the arguments u, k. 

Then, from the integral equations for and Wi^ we notice the relationship 

ik dv h2 . . 

which we use to eliminate Wi t 2 from these equations; so, we are left with the 
second set of two integral equations in t>i j2 : 

for z > d/2 

C 2 K 2_ U 2 vi-- ^a^) J d/2 az e 'v^z )- 

i^l^-^l) r-d/2 j , i K f z - z ') / a, / 9 o\ 

+ 2(c 2 K 2 -^) e v^a/z) 2 ( c 2 K 2_ w 2)e 6> 2 ^ a/z; 

and 

-iS^rf 2 ^ 1 -' 1 ^^)- (24) 

- 27 S^ ye -^- d /% 1 (rf/2) + ^g^e--(^/%2(-rf/2) 

for 2! < —d/2. It is worth observing in deriving these equations the non- 
intervertibility of the derivatives and the integrals, according to the identity 

£ / dz'f(z)£re lK \ z - z I = / ^7(^> m l Z ~ 2 ' - 2inf(z) (25) 

°^ Jd/2 °^ Jd/2 
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for any function f(z), z > d/2; a similar identity holds for z,z' < —d/2. 

It is due to the discontinuity in the derivative of the function eH z z I for 
z = z . We can see that these equations become equations (JTOj) in the non- 
retarded limit by taking formally the limit c — > oo. However, this is not so 
for their dispersion equations, as we shall see below. One can also see from 
equations (1211) . (1231) and (1241) that the coupling between the two bodies is 
performed through both bulk and surface degrees of freedom, in contrast to 
the non-retarded situation, where this coupling occurs only through surfaces 
(equations (ITT]) ). 



We turn now to equations (T2TT) . Taking the second derivative with respect 
to z in these equations we get 



92v h2 , ( .2 ^1,2 



+ U 2 - -M U-l = o , (26) 



dz* T \ c 2 ; ^ 

which tells that v± 2 are a superposition of two waves e ±lKl • 2Z , where 



« 1)2 = xh 2 -'-^ ■ ( 27 ) 



c 2 



We note that such modes are polaritonic modes, since u 2 = c 2 (k 2 + k 2 ) = 
c 2 (k 2 + k 2 2 ) +w 2 2 = c 2 K 2 2 + cuf 2 , where K 1)2 = (k, k 1i2 ), which is the well- 
kown dispersion relation for the polaritonic modes. It can also be written 
as cu 2 £i,2 = c 2 Kf 2 , where £1,2 = 1 — oj\ 2 /uj 2 is the dielectric function for 
metals. This relation is well-known in the so-called thery of "effective medium 
permittivity". We take v± 2 = A^ 2 e lKiaz ', where A 1>2 are amplitudes to be 
determined. Then, equations (I2TI) have non-trivial solutions for frequencies 
u given by the roots of the dispersion equation 



w («i + k)(k 2 -«) 



(«i - k)(k 2 + k) 



(28) 



Equation (1251) has a branch of roots for the damped regime (evanescent 
modes) Ki = ia\, k 2 = — iac2, given by 

tan^= " ( 2 ai + a2) , (29) 

where 



«1,2 = \/ - K 2 , ^1,2 > , (30) 
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and k real. Since these modes are damped inside the bodies and propagating 
in-between the bodies they may be called surface plasmon-polariton modes. 
It is worth noting the correct choice of the sign of the square root in this case, 
in order to get the correct behaviour at infinity, = A\e~ axz for z > d/2 
and v 2 = A^e" 22 for z < —d/2. The roots of equation (j29i) can be written 
as 

tt 2 x 2 



^ = 0^0 + -^, (31) 

where x = and n — 1/2 < x n < n + 1/2, n = 1,2,3,... for x n < 
min (uji, u 2 ) d/irc. For identical bodies the roots are given by 



7T 2 T) 2 

n = cyjk* + — (32) 

for any integer n = 0,1,2.... They correspond to propagating (polariton) 
modes (ki = n 2 and k all real numbers) and arise from equation (1281) for 
e 2ind _ 1 Equation fl29l) may have another solution in the vicinity of the 
vertical asymptote of the function in its rhs, which, however, is irrelevant for 
our discussion. 



Similarly, 1*1,2 from equations (1231) and (l24|) obey the same equation (|26l) . 
We look again for solutions of the form v 1>2 = A^ 2 e lKl ' 2Z , where A 1>2 are 
amplitudes to be determined. According to equations (122]) these modes are 
transverse modes, as they should be (for Ki )2 real). The relevant dispersion 
equation is given by 

e 2iKd = K + K,)(k2 - + k 2 ){KK 2 - k 2 ) 

(ki — k)(k 2 + k)(kki — k 2 ){KK 2 + k 2 ) 

We note that this dispersion equation does not become the non-retarded 
dispersion equation (1281) by taking formally the limit c — > 00. 



An analysis similar to the one performed above for equation (12811 shows that 
equation (1331) has a branch of roots 



n 2 = c^k 2 + ^jf, (34) 

where y = and y n < min (u>x, uj 2 ) d/nc. They correspond to surface 
plasmon-polariton modes Ki = icti,K 2 = —ia 2 and k real. We note that 
y n may differ from x n . For identical bodies these roots are those given by 
equation (132|) . Some other isolated roots may appear, as for instance the 
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one corresponding to an overall damping, i.e. K\ = iai, k 2 = — 2a; 2 , K = ia, 
where a = \Jk 2 — lu 2 /c 2 , uj < ck. It is given by 



where min (u>i, u 2 ) < ir^/2cz /d < max (ui, uj 2 ). Such an isolated mode does 
not contribute significantly to the energy, so we may neglect it in our subse- 
quent analysis. 

We can take the limit d — > oo in equation (1331) . It can be shown that this 
limit amounts formally to put e 2%Kd = 0. [23] We get in this case the surface 
plasmon-polariton modes corresponding to a semi-infinite body, given by 
aaci 2 = k 2 , i.e. 

2uj 2 c 2 k 2 

uj 2 = 1,2 , (36) 

ul 2 + 2c 2 k 2 + Ju\ 2 + 4c 4 fc 4 

as derived previously. [2^ In general, there are problems with taking formally 
the limits d — > or d — > oo in the above equations, as expected. 

It is also worth interesting to look for solutions of the type 

vi, 2 = A h2 [e iKl ' 2Z - e ±iK ^^ z) ] (37) 

for equations ([231 and (l24]h which are vanishing on the surfaces, v\p{±d/2) 
("fixed surfaces" boundary conditions). In this case, we get again the res- 
onance modes Q given by equation (1321) . irrespective of the bodies being 
distinct or identical. In addition, we may get special modes u = o;i j2 , 
uj 1 = c 2 k 2 + uj\ 2 (aci ( 2 = 0) or uj = ck (k = 0), which do not depend on 
distance d. Other boundary conditions can be put on surfaces z = ±d/2, 
and we can get the corresponding eigenmodes. 

We note that the dispersion equations (1281) and (1331) appear, though in a dis- 
guised form, in various formulations of the fluctuations theory. [2 j . [3 J . [5] . [32] 
Within the framework of this theory the dielectric function is included from 
the beginning. On the contrary, we recover the dielectric function in the final 
results of the present approach, which shows that our approach is equivalent 
with the so-called "effective medium permittivity" theory. 

We pass now to the zero-point energy corresponding to the f2i 2 -eigenmodes 
given by equations (13T1) and or the f2-branch given by equation (1321) 

(for identical bodies or "fixed surfaces"), in the limit min (u>i, u; 2 ) d/itc 3> 1. 
These are the only eigenfrequencies which depend on distance d. In the limit 
min (cg>i,cg> 2 ) djuc 1 these modes are dense sets, and it is easy to see that 
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their contributions to the zero-point energy are equal (corresponding to the 
two polarizations), so we can write the total zero-point energy as 



/ rT7-2 /">"' 2 

E = fc Ef 2 + ir. (38) 

kn=0 

where defined above; for identical bodies (or for " fixed surfaces") 

x n = n. We follow the standard regularization procedure by removing the 
ultraviolet divergencies and using the Euler-MacLaurin formula. |34| As it is 
well-known, the energy thus regularized reads 

where B 2 k are Bernoulli's numbers and 



f(x)= I dk.k\lk* + 7 ^= l - [ du-^x. (40) 

JO V (1 Z J n 2 x 2/ d 2 

Since x = (and y — 0), we get the well-known energy E = —^hcB^/Ald 3 = 
—ii 2 hc/720d 3 and Casimir force F = 7r 2 ftc/240<i 4 per unit area. The same re- 
sult is obtained for the f2-modes given by equation (1321) with n = 0, 1, 2..., cor- 
responding to identical bodies or the "fixed surfaces" boundary conditions 
t>i ;2 (±<i/2). It is easy to see that for decreasing min (o>i, u^) d/irc the number 
of x„-roots contributing to energy decreases, the numerical coefficient of the 
Casimir force decreases gradually, and the <i~ 4 -dependence deteriorates, untill 
a cross-over may occur to the non-retarded van der Waals-London d~ 3 -force. 



The dispersion equations (12811 and (1331) hold also for dielectrics, providing 
the wavevectors K\p are changed according to 



.2 .2 W 1.2 



<2 - *' - — • (41) 



c u> — uj 012 



We can get a usual model of dielectric for cjoi,2 ^ ^1,2 • In this case, the 
wavevectors k± 2 become 



uj{ 2 UJ 2 



^2 = \h 2 + -T ± -w , (42) 



wgi 2 c 2 



and we cannot have anymore surface plasmon-polariton modes (evanescent 
modes). In general, under these circumstances, the dispersion equations (1281) 
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and (1331) have no solutions, except for identical bodies when we may have the 
f2-modes given by equation f|32l) (e 2iKd = 1) for n = 0,1,2.... These modes 
correspond to propagating polaritons and give again the classical result for 
the Casimir force F = 7r 2 hc/240d A per unit area. Similarly, for a dielectric- 
metal pair there is no force, except for boundary conditions when 
the resonant f2-modes given by equation (1321) for n = 0,1,2... are present. 
The latter result holds for any pair of bodies. It is, however, worth stressing 
that such results depend on our model of dielectric function for dielectrics, 
and, in general, it is necessary to have a quantum-mechanical treatment for 
the internal dynamics of the dielectrics. 



4 Discussion and conclusions 



In conclusion, we may say that we have derived here van der Waals-London 
and Casimir forces acting between two semi-infinite bodies with parallel sur- 
faces by calculating the electromagnetic eigenmodes in matter and estimating 
their zero-point energy (vacuum fluctuations). We have adopted well-known, 
simple, usual models for matter polarization in metals and dielectrics and 
made use of the equation of motion for the polarization in order to get cou- 
pled integral equations. The eigenfrequencies of these equations have been 
identified and used in calculating the zero-point energy. In the non-retarded 
(Coulomb) limit we get the well-known van der Waals-London <i~ 3 -force, 
arising from the surface plasmons, where d is the distance between the two 
bodies. The numerical coefficient of this force acquires various values, de- 
pending on the nature of the bodies and on their being distinct or identical. 
When retardation is included we get the Casimir <i~ 4 -force arising from sur- 
face plasmon-polariton modes (evanescent modes) for a pair of metals. The 
classical numerical coefficient of this force (7r 2 /240) is obtained for distances 
much larger than the characteristic wavelengths (~ c/ooi^, where are the 
plasmon frequencies) of the bodies, and it diminishes gradually for shorter 
distances, while the force loses its characteristic d~ 4 -dependence. For a pair 
of identical dielectrics we get the classical Casimir result arising from prop- 
agating polariton modes. The same result holds for any pair of bodies with 
"fixed surfaces" boundary conditions. 

As it is well-known, the fluctuations theory[32] predicts Casimir forces be- 
tween any pair of bodies, in contrast with our results, which give a vanishing 
force for two distinct dielectrics, for instance. The difference originates in 
the circumstance, usually overlooked, that the equivalent of our dispersion 
equations (1281) and (133]) in the fluctuations theory have no solutions in some 
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for instance, for distinct dielectrics. The usual theorem of meromor- 
phic functions, applied within the framework of the fluctuations theory, [4j-[6j 
gives then a finite result, but it does not represent the energy of the eigen- 
modes. The problem does not appear in the non-retarded regime, where our 
results coincide with those of the fluctuations theory. On the other hand, 
we must stress again upon the fact that our model for the dielectric func- 
tion may not be perfectly adequate for describing the internal polarization 
of dielectric matter. Again, this is immaterial in the non-retarded regime, 
and we succeeded to compute a <i~ 4 -van der Waals-London force between 
a classical model of polarizable point-like particle and a semi-infinite body. 
But our approach fails in this case in the retarded regime, where a quantum 
mechanical treatment is necessary, as in the original attempt in Ref. [35]). 

Finally, it is worth noting that the dispersion equations (1251) and (1331) can also 
be obtained by calculating the reflected field in-between the bodies (fields for 
semi-infinite bodies). [23] If fi,2 are the amplitudes of these fields (for a given 
polarization), then the dispersion equations (1251) and (1331) are obtained from 
ri = r 2 e 2lKd . We note that |ri )2 1 2 are the reflection coefficients, and for two 
perfectly reflecting bodies \ri\ — |r 2 | = 1. If we neglect the phases of the 
coefficients r 1>2 , and put r\ — r 2 = 1, we get the Casimir dispersion equation 
e 2iKd _ ^ (Q-modes given by equation (132|) ). However, it is precisely these 
phases that give the damped surface plasmon-polariton regime, as we have 
shown in the present paper, and these phases are not equal in the damped 
regime, not even for identical bodies. This is related to the correct choice 
of the sign of the square root in Ki i2 , which, as we have shown here, is 
Ki = iai and k 2 = — ia 2 (equations (1291) and (1301) ). For the propagating 
regime (vanishing phases) and identical bodies (ri = r 2 ) we get again the 
Casimir dispersion equation e 2tftd = 1, as we do for "fixed surfaces" boundary 
conditions (in the latter case irrespective of the bodies). 
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